On the energy-momentum tensor for a scalar field on manifolds 

with boundaries 



Aram A. Saharian* 
Department of Physics, Yerevan State University, 1 Alex Manoogian Str., 

375049 Yerevan, Armenia, 
and 

the Abdus Salam International Centre for Theoretical Physics, 

34014 Trieste, Italy 

February 1, 2008 



Abstract 

We argue that already at classical level the energy-momentum tensor for a scalar field on 
manifolds with boundaries in addition to the bulk part contains a contribution located on 
the boundary. Using the standard variational procedure for the action with the boundary 
term, the expression for the surface energy-momentum tensor is derived for arbitrary bulk 
and boundary geometries. Integral conservation laws are investigated. The corresponding 
conserved charges are constructed and their relation to the proper densities is discussed. 
Further we study the vacuum expectation value of the energy-momentum tensor in the cor- 
responding quantum field theory. It is shown that the surface term in the energy-momentum 
tensor is essential to obtain the equality between the vacuum energy, evaluated as the sum 
of the zero-point energies for each normal mode of frequency, and the energy derived by the 
integration of the corresponding vacuum energy density. As an application, by using the 
zeta function technique, we evaluate the surface energy for a quantum scalar field confined 
inside a spherical shell. 



PACS numbers: 03.50.Kk, 03.70.-Fk, 04.62. -hv, 11.30.-j 



1 Introduction 

In many problems we need to consider the physical model on background of manifolds with 
boundaries on which the dynamical variables satisfy some prescribed boundary conditions. In 
presence of boundaries new degrees of freedom appear, which can essentially influence the dy- 
namics of the model. The incomplete list of branches where the boundary effects play an 
important role includes: surface and finite-size effects in condensed matter theory and statisti- 
cal physics, theory of phase transitions and critical phenomena [1, 2], canonical formulation of 
general relativity (see [3]), the definition of the gravitational action, Hamiltonian and energy- 
momentum [4]-[16], path-integral approach to quantum gravity and the problem of boundary 
conditions for the quantum state of the universe (see, for instance, [17] and references therein), 
bag models of hadrons in QCD, string and M-theories, various braneworld models, manifolds 
with horizons and the thermodynamics of black holes, AdS/CFT correspondence [18], holo- 
graphic principle [19, 20]. In recent years, the study of the boundary effects in quantum theory 
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has produced several important results. For example, an interesting idea is to consider the black 
hole event horizon as a physical boundary [21]. This induce an extra term in the action, having 

as consequence the existence of a central charge in the algebra of generators of gauge transfor- 
mations. Using this central charge it is possible to determine the asymptotic behavior of the 
densities of states and in this way to get the entropy for a black hole [22]. In the context of 
string theory, the D-branes are natural boundaries for the open strings, having very interesting 
effects in the theory. 

In quantum field theory the influence of boundaries on the vacuum state of a quantized field 
leads to interesting physical consequences. The imposition of boundary conditions leads to the 
modification of the zero-point fluctuations spectrum and results in the shifts in the vacuum 
expectation values of physical quantities, such as the energy density and vacuum stresses. In 
particular, vacuum forces arise acting on constraining boundaries. This is the familiar Casimir 
effect. The particular features of the resulting vacuum forces depend on the nature of the quan- 
tum field, the type of spacetime manifold and its dimensionality, on the boundary geometries 
and the specific boundary conditions imposed on the field. Since the original work by Casimir in 
1948 [23] many theoretical and experimental works have been done on this problem, including 
various types of bulk and boundary geometries (see, e.g., [24]-[31] and references therein). The 
Casimir force has recently been measured with a largely improved precision [32] (for a review 
see [28, 33]) which allows for an accurate comparision between the experimental results and the- 
oretical predictions. An essential point in the investigations of the Casimir effect is the relation 
between the mode sum energy, evaluated as the sum of the zero-point energies for each normal 
mode of frequency, and the volume integral of the renormalized energy density. For scalar fields 
with general curvature coupling in Ref. [34] it has been shown that in the discussion of this ques- 
tion for the Robin parallel plates geometry it is necessary to include in the energy a surface term 
concentrated on the boundary (see Refs. [35, 36] for similar issues in spherical and cylindrical 
boundary geometries and the discussion in Ref. [37]). However, in Refs. [34, 35, 36] the surface 
energy density is considered only for the case of the Minkowski bulk with static boundaries. In 
the present paper, by using the standard variational procedure, we derive an expression of the 
surface energy-momentum tensor for a scalar field with general curvature coupling parameter 
in the general case of bulk and boundary geometries. It is well-known that the imposition of 
boundary conditions leads to a vacuum expectation value of the energy-momentum tensor with 
non-integrable divergences as the boundary is approached [38, 39, 40]. The role of surface ac- 
tion to cancel these divergencies is emphasized in Ref. [40], where it has been shown that the 
finiteness of the total energy of a quantized field in a bounded region is a direct consequence of 
renormalization of bare surface actions. How boundary corrections arc implemented by surface 
interaction and the structure of the corresponding counterterms in the renormalization proce- 
dure are discussed in Ref. [41] within the Schrodinger representation of quantum field theory. 
By suitably modifying the classical action with boundary corrections, in Ref. [42] it has been 
shown that for scalar fields confined in a cavity familiar functional methods can be applied and 
the functional measure is constructed by treating the boundary conditions as field equations 
coming from a variational procedure. Surface terms play an important role in the recently pro- 
posed procedure [11] (see also [12]-[15] and references therein) to regularize the divergencies of 
the gravitational action on noncompact space. By adding suitable boundary counterterms to 
the gravitational action, one can obtain a well-defined boundary energy-momentum tensor and 
a finite Euclidean action for the black hole spacetimes. 

This paper is organized as follows. In Sec. 2 the notations are introduced for the geometrical 
quantities describing the manifold and the structure of the action is discussed. In Sec. 3 the 
variations of the bulk and surface actions are considered with respect to variations of the metric 
tensor and the expressions for the bulk and surface energy-momentum tensors are derived. The 
integral conservation laws described by these tensors are studied in Sec. 4. The expressions 



2 



for the corresponding energy-momenta are derived and relations to the proper densities are 
discussed. The vacuum expectation values of the energy-momentum tensor for a scalar field on 
a static manifold with boundary are considered in Sec. 5. It is argued that the surface part 
of the energy density is essential to obtain the equality between the vacuum energy, evaluated 
as the sum of the zero-point energies for each normal mode, and the energy obtained by the 
integration of the corresponding vacuum energy density. As an application of general foemula, 
in Sec. 6 we consider the surface energy for a quantum scalar field confined inside a spherical 
shell. Section 7 concludes the main results of the paper. An integral representation of the partial 
zeta function needed for the calculation of the surface energy on a spherical shell is derived in 
Appendix A. 

2 Notations and the action 

Consider a D + 1 - dimensional spacetime region M with metric go. and boundary dM. We 
assume that the spacetime M is foliated into a family of spacelike hypersurfaces S defined by 
t = const. The boundary of M consists of the initial and final spacelike hypersurfaces Si and 
E2, as well as a timelike smooth boundary dMg and, hence, dM = SiUS2U3Ms. The boundary 
of each S we will denote dT,: dT, = S n dMg. The corresponding spacetime region is depicted 
in Fig. 1 for a three dimensional M. The inward pointing unit normal vector field on dM is 
denoted n*. This vector is normalized in accordance with rijn* = e, where e = 1 and —1 for 
spacelike and timelike boundary elements, respectively. We define n* as the future pointing unit 
normal vector to the spacelike hypersurfaces S. On the hypersurfaces Si and S2 this vector is 
connected to the vector by relations: = on Si and u' = — on S2. In the consideration 
below we will not assume that the timelike boundary dMg is orthogonal to the foliation of the 
spacetime. We define as the forward pointing timelike unit normal vector field to the surfaces 
dT, in the hypersurface dMg. By construction one has u'^rii = 0. We further define n* as the 
vector field defined on dMg such that is the unit normal vector to dT, in hypersurface S and, 
hence, by construction n*nj = 0. It can be seen that the following relations take place between 
the introduced vectors (see, for instance, [8, 10, 16]) 

h' = X {n' - rju') , u' = X {u' + rin') , (2.1) 

where is understood as the normal to dMs, the variable 77 = u^rii measures the non- 
orthogonality and A = 1/ \/l + rf . For rj = the foliation surfaces are orthogonal to the 
boundary dMg. Note that the quantity A is connected to the proper radial velocity v by the 
relation [16] A = Vl — f^. 

In the present paper we deal with a real scalar field ip{x) propagating on the manifold M. 
The action describing the dynamics of the field we will take in the form 

S = Sb + Ss, (2.2) 

where the first term on the right, 

Sb = ^ [ d^'+'x^l (V.^'VV - - CRV'^) , (2.3) 

is the standard bulk action for the scalar field with curvature coupling parameter R is the 
Ricci scalar for the manifold M, is the covariant derivative operator associated with the bulk 
metric gn^ (we adopt the conventions of Ref. [43] for the metric signature and the curvature 
tensor). The second term on the right of Eq. (2.2), 

Ss = -e f d'^x^/\h\ {((p^K + ms(p^) , (2.4) 
JdM 
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Figure 1: Manifold M with boundary dM = Ei U E2 U dMg fohated by spacehke hypersurfaces 
E with boundaries dT,. n* and are unit normal vector fields to the hypersurfaces dMg and E 
respectively, is the unit normal vector to the {D — l)-surface dH in the hypersurface S, and 
is the unit normal vector to 5E in the hypersurface dMg. 

is the surface action with a parameter m^. In Eq. (2.4), h is the determinant of the induced 
metric hik for the boundary dM, 

hik = Qik - euiUk, (2.5) 

which acts as a projection tensor onto dM, K is the trace of the extrinsic curvature tensor Kik 
of the boundary dM: 

K = g'''Kik = h'''Kik, Kik = hlhl'Vinm, h\ = g''hik. (2.6) 

For the most important special cases of minimally and conformally coupled scalars in Eq. (2.3) 
we have ( = 0,{D — 1)/4D, respectively. Note that in Eq. (2.4) the term with the extrin- 
sic curvature scalar corresponds to the Gibbons-Hawking surface term in general relativity, 
— e Jgf^ d^x^/\h\ K. This term is required so that the action yields the correct equations of mo- 
tion subject only to the condition that the induced metric on the boundary is held fixed. Note 
that for non-smooth boundary hypersurfaces it is necessary to add to action (2.2) the corner 
terms (for these terms in the gravitational action of general relativity see [7, 9, 10, 16]) with 
integrations over {D — 1) -hypersurfaces at which the unit normal changes direction discon- 
tinously. As we have assumed a smooth boundary dMg, in our case these hypersurfaces are 
presented by the {D — l)-surfaces 9Ei = Ei n dMg and 9E2 = E2 fl dMg and the corresponding 
corner term in the action has the form 

Sc = -C / d^-^x^\(p'^avcsiiLh{r)), (2.7) 

where Jg^^ = Jgj.^ — Jgj.^ and a is the determinant of the metric tensor 

= gik + nink - UiUk = gik + nifik - UiUk (2.8) 

induced on dH. In this paper we are interested in the variation of the action with respect to 
the variations for which the metric and field are held fixed on the initial and final spacelike 

hypersurfaces Si and E2, and consequently on OSi and 9E2. In this case the corner terms 
will not contribute to the variations under consideration and we will omit them in the following 
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consideration (see, however, the discussion after formula (3.27) in Sec. 3 for the total divergence 
term in the variation of the action, located on the boundary). Note that the contribution of the 
corner terms will be important in the case of a non-smooth boundary dMg having wedge-shaped 
parts. 

3 Variation of the action and the energy-momentum tensor 

3.1 Variation of the bulk action and the volume energy- momentum tensor 

In this subsection we consider the variation of the bulk action (2.3) induced by an infinitesimal 
variation Sgik in the metric tensor. As the careful analysis of the surface terms is crucial for our 
discussion we will present the calculations in detail (for the corresponding calculations in general 
relativity see, for instance, [16]). To find the expression for the energy-momentum tensor of the 
scalar field with action (2.2), we need the variation of this action induced by an infinitesimal 
variation Sgik in the metric tensor. 

First of all we consider the variation of the term involving the Ricci scalar: 

Sg f d^+'x^\ip^R= f d''+^x^/\i\vH^9"'Gik + g"'SRik), (3.1) 

JM JM ^ ' 

where and Gj^ = — gi^^Rjl arc Ricci and Einstein tensors, respectively. To evaluate the 
second term in the braces on the right we note that the contracted variation g^'^SRik is a pure 
spacetime divergence, 

g"'6Rik = -ViW\ W'=g"'6T\k-g''6Ti, (3.2) 

with r'j^. being the Christoffel symbol. In order to find the variations of the Christoffel symbol 
in this formula, note that from the identity S{Vigki) = it follows 

g^'Srl = -gHri + g^'^g'PVpSgik. (3.3) 

Combining this with the relations 

dVi = ViiSln ^\) = -^gi^ViiSg'n, q'^ ff^^ p^9ik = -V^(V"*), (3.4) 

one finds 

W = ff''fffc™V,(55'") - V„,( V)- (3.5) 
By taking into account formula (3.2) and using the Stoke's theorem 

/ d''+'x^M^^V' = -e f d''x^/\h\n^V\ (3.6) 
Jm JdM 

for the term with the variation of the Ricci tensor one has 



/ d°+^x^/\^\^'^g'^5Rik= f d^+^xy/\i\W^Vi^^ + e [ d^xVW/n^W'. (3.7) 
Jm Jm JdM 

For the first term on the right, using formula (3.5), we find 

/ d^'+'xVWlW'Vi'p^ = [ d''+'x^\dg'''(ViVk'p'-g'"'gikViVm^^) 

JM JM ^ ' 

+e / d''x^\Sci^ (niVfe¥^2 _ g.^n^^^^A . (3.8) 

JdM ^ ' 
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Substituting this into formula (3.7) and taking into account formula (3.5), for the term containing 
the variation of the Ricci scalar, Eq. (3.1), we obtain the formula 



JM JM ^ ' 



JdM ^ ^ 

'9kmn'Vi{Sg''^)-niVm{Sg'^)]}. 
This leads to the following variation of the bulk action 

S,s, = U d^^'xVW\Sg''Tt'^ + U d^xVWlUg'H-CriiVw' 

^ JM ^ JdM 

'nlVUSg'"')-9kmn'^^{6g'n]}, 



(3.9) 



(3.10) 



where is the standard metric energy-momentum tensor for a scalar field with a curvature 

coupling parameter ( (see, for instance, [43]): 



.(vol) 



m 



-ik - ^i^^k^-2^ikg^^l'P^m<P+^<p'9ik-C^'Gik + C9ik9''^'^Nm^-C^iVkV^ (3.11) 
Note that this expression can be also written in the form 

where 



(3.12) 



'-ik 



^iV'^kV + 



C - ^ ) gikg^'^yNm - CRik - CV^V^ 



(3.13) 



This tensor does not explicitly depend on the corresponding mass. For the covariant divergences 
of tensors (3.11) and (3.13) it can be seen that 



(vol)i 



^ (^''"ViV^v' + CR^ + "iV) 



(3.14) 
(3.15) 



From formula (3.10) it follows that the surface term in the variation of the bulk action contains 
the derivatives of the variation for the metric. In next subsection we will see that these terms 
are cancelled by the terms coming from the variation of surface action (2.4). 



3.2 Veiriation of the surface action and the surface energy-momentum tensor 

Now we turn to the evaluation of the variation for surface action (2.4) induced by the variation 
of the metric tensor. This variation can be written in the form 



SgSs 



-e / d^x^/\h\ 

'OM 



First of all we note that for the variation of the extrinsic curvature tensor trace one has 

5gK = dgig'^^ViUk) = Sg'^ViUk+g'^'Vidnk - nig'HT\k. 



(3.16) 



(3.17) 
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By making use formulae 



Sui = ^enin''n''dgki, ViUk = Kik + enin'''Vink, 



(3.18) 



and formula (3.3) for the contracted variation of the Christoffel symbol, it can be seen that 



6gK = -K^5gik + ^Kn'n^5gik + n'hl'^ ( \Vi5gki - VkSgu 



(3.19) 



2 \2 

Substituting this relation into formula (3.16), for the variation of the surface action one finds 



SgSs 



JdM 



(3.20) 



(3.21) 



By taking into account that in Eq. (3.10) we can replace 

niVmiSg'"^) - gkmn'Vi{5g^n = g^^^ri' NiiSgkm) - ^^miSgik)] , 
and using formulae (3.10) and (3.20), for the variation of the total action one finds 

^ Jm ^ JdM 

+ Cgiku'Viif^ - (ip^ {2K,k - Khik - en^rikK)] + C^'n^/i'" (Sgi^)] . (3.22) 

Let us consider the last term in the figure braces on the right which is the only term containing 
the derivatives of the metric tensor variation on the boundary. To transform this term, we write 
it in the form 

n'h'"'Vl{Sgim) = '^l{n'h'"'Sgim) - Sgim'^iin'h'^). (3.23) 
To proceed further we note that for an arbitrary vector field V'' one has 



V,y' = DiV^ + en'niViV^ = DiV^ + en'Vi{niV^) - en'V^ViUi, 



(3.24) 



with Di being the covariant derivative projected onto dM. This allows us to write the first term 
on the right of Eq. (3.23) in the form 



Viin'h'^Sgim) = Diin'h'^Sgim) - en'n'h'^5gkrrN m, 



•Im 



(3.25) 



where we have taken into account that nih''"^ = 0. Now using the second formula in Eq. (3.18), 
one finds 

n'h^'^Viidgim) = Di{v}h}^5gim) - Sgik {k'^ - en'n'^K^ . (3.26) 

Substituting this into formula (3.22), we receive to the final expression for the variation of the 
total action 

SgS = I [ d''+'x^MSg"'T^r'^ + ^ f d''x^/\h\\Sg"'[hikmsv'' 

^ JM ^ JdM ^ 

+ ChikU^Vi^^ - C^' {Kik - Khik)] + CA(^'"*/i'"'^5im)} • (3.27) 

An important feature of this form is that the derivatives of the metric tensor variation are 

contained in the form of a total divergence on the boundary (the last total divergence term in 
the figure braces). To see the relation between this total divergence term and the variation of 



the corner term (2.7) in the action, note that using formula (3.18) for the variation Srii, it can 
be seen that hJ-'^n^dgim = —h'j^Sn'^. With this replacement, decomposing the integral with the 
total divergence term into integrals over dMg, Si, S2 and using the Stoke's theorem, we can see 
that 

^ JdM ^ J(9Ei 

d^-^a:i/|CT|(^25arcsinh(r?). (3.28) 

In obtaining the second equality we have used relations (2.1) (for similar calculations in general 
relativity see, for instance, [16]). Now we see that the total divergence term is cancelled by the 
corresponding term coming from the variation of the corner part of the action, Eq. (2.7). 

Now we return to the variation given by Eq. (3.27). Using the Stoke's theorem, we can 
transform the integral with the total divergence term into the integrals over {D — l)-surfaces 
9Ei and We consider variations for which the metric is held fixed on the initial and final 
spacelike hypersurfaces Si and E2 and, hence, on dYl\ and dYl^- Consequently, for the class of 
metric variations Sg^^ with (55**^|s, = 0, j = 1,2, the variation of the action takes the form 



\ [ d^+'x^Wl^g^'Tt,"'^ + \l d^x^ , (3.29) 



where we have introduced a notation 

Tik = C^'^Kik-hik(Q'/K + (:n^Vl^'^ + ms^'^') , (3.30) 

and have taken into account that e = — 1 for the timelike hypersurface dMg. 

Using the definition of the metric energy-momentum tensor as the functional derivative of 
the action with respect to the metric, one finds 



_ ^(vol) ^ ^(surf) 



J-ik - /r-C X^ik ~ -^ik '^-^ik ■ I^O.OIJ 



Here the part of the energy-momentum tensor located on the surface is defined by relation 

tI^''^ =5{x;dMs)T,k, (3.32) 
and the 'one-sided' 5-function 5{x;dMs) is defined as 

f d^^^x^\5{x-dMs)= ! d^x^/\h\. (3.33) 

An important property of the surface energy-momentum tensor is that it is orthogonal to the 
normal vector field n*: 

which directly follows from the expression for the tensor Tik- Another important point is that, 
the surface energy-momentum tensor does not depend on the foliation of spacetime M into 
spacelike hypersurfaces S. Using expressions (3.11) and (3.32) for the volume and surface 
energy-momentum tensors and the Gauss-Codacci equation 



h^ViK = hl^h'^'ViKik + hl^n'Rik, (3.35) 
the following relation can be derived 

Dkri - Ukh^^T^^'''^'' = (4V;rf - n^T^^™')*^) = -h^ {2C^K + 2ms^ + n'Viip) V^^ (3.36) 
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on the boundary dMg. 

Note that in all relations obtained so far we have not used the field equation for ip{x). Now 
we turn to this equation and corresponding boundary condition for the field following from the 
extremality of the action with respect to variations in the scalar field. The variation of action 
(2.2) induced by an infinitesimal variation Stp in the field has the form 



S,pS 



Jm 

-e d'^x^/\h\^(p(2((pK + n^Vi(p + 2ms(p) . 

JdM ^ ' 



(3.37) 



For the variations with = on the initial and final hyper surf aces Si and S2, this leads to 
the following field equation on the bulk 



and to the boundary condition of Robin type, 

2(Cif + ms)(^ + nV;<^ = 0, 



(3.38) 



(3.39) 



on the hypersurface dMg- Note that the Dirichlet boundary condition is obtained from here in 



the limit — > 00. 



Prom relation (3.12) for the tensors T^"^^ and T\^^ we see that on the class of solutions to 



(1) 



(vol) 



T^A^ . An advantage of the form TV"* in 



(1) 



field equation (3.38) these tensors coincide: 7^. 
calculations of the vacuum expectation values of the energy density in the corresponding quan- 
tum theory is that it is sufficient to evaluate only two quantities, (0|(/?^(a;)|0) and (0|(9o(/?)^|0). 
In the case of T^"^ ■, due to the second term on the right of Eq. (3.11), we need the vacuum 
expectation values for all derivatives (0|5;<p5^<^|0). In fact, this is the reason for our choice of 
the form in investigations of the Casimir densities for various types of geometries. Note 
that, using boundary condition (3.39), the surface energy-momentum tensor can be also written 
in the form 



p(surf) 



6{x; dMs 



ik 



2(-'- 

^ 2 



hik(pn''Vnp 



(3.40) 



The surface energy-momentum tensor in this form does not explicitly depend on the parameter 
nis. From relations (3.14) and (3.36), by making use field equation (3.38) and boundary condition 
(3.39), we see that on the solutions of the field equation the volume energy-momentum tensor 
covariantly conserves. 



(vol)i 



0, 



and for the surface energy-momentum tensor one has 



(3.41) 



(3.42) 



The last equation expresses the local conservation of the boundary energy-momentum tensor up 
to the fiow of the energy-momentum across the boundary into M. Integral consequences from 
formulae (3.41) and (3.42) we will consider in Sec. 4. 

The traces for the volume and surface energy-momentum tensors are easily obtained from 
(3.13) and (3.40): 



(3.43) 

(3.44) 
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where Q = (D — l)/4-D is the curvature coupUng parameter for a conformally coupled scalar 
field. Hence, for a conformally coupled massless scalar with nis = both these traces vanish. 
The latter property also directly follows from the conformal properties of the action (2.2). 
Indeed, it can be seen that under conformal transformations of the metric gn^ = il'^{x)gik and 



the corresponding transformation of the field f 
transforms as 



0(i-'D)/2y,^ action (2.2) with m 



0, 



2</?VVVi Inf^ - 



D - 1 



S[gik, (p] = S[gik, ^\ + D{Q - Co) / d''+'x^\g\ 

JM I ^ 

(3.45) 

and, hence, is conformally invariant in the case ^ = Cc- In order to obtain (3.45) we have used 
the standard transformation formula for the Ricci scalar and the corresponding formula for the 
extrinsic curvature tensor: 

= ^K,^ + h^kU^VSi . (3.46) 

From the conformal invariance of the action it follows that boundary condition (3.39) is also 
conformally invariant when nis = and C = Cc- In this case it coincides with the conformally 
invariant Hawking boundary condition [40] . This is also directly seen from the relation 



2C,K(p + n'ViCp = n 



D+l 
2 



2CKip + nWiip + 2D{C - Qfn'Vi InO 



(3.47) 



where = ^n} . Note that for general m^, under the conformal transformations the surface 
energy-momentum tensor transforms as 

fl''''''^\^{x)] = S{x;dMs)fn^{x)] 

= n-''-'{T^'''''^''[^{x)]-D{C-Cc)S{x;dMs)h^^''n'Vi\nn], (3.48) 

where S{x; dMg) = Q,~^6{x; dMg). For C = Cc this transformation formula is the same as for the 
volume part of a conformally coupled massless scalar field. 

As an example let us consider the special case of a static manifold M with static boundary 
dMs. In an appropriately chosen coordinate frame one has gop = and the components goo, 
gfSp, P,p = 1,2,...,D, do not depend on the time coordinate. Taking into account that in 
static coordinates nP = 0, for the 00-component of the extrinsic curvature tensor we find Koo = 
^nl^d/sgoo- Substituting this into Eq. (3.40), for the surface energy density we obtain the 
expression 



fp(surf)0 



S{x;dMs) 



((f'^n'^ d/siln y/g^) 



Q^H{x- dMs) l^Cnf^df^iln^) + (4C - l)(i^^ + m^/C) 



(3.49) 



where in the second line wc have used boundary condition (3.39). In the particular case goo = 1 
the first expression coincides with the formula of the surface energy density used in Ref. [40] 
to cancel the divergence in the vacuum expectation value of the local energy density for a 
quantized scalar field and in Refs. [34, 35, 36] to evaluate the vacuum expectation values for 
flat, spherical and cylindrical static boundaries on the Minkowski background. Note that in 
this case the surface energy vanishes for Dirichlet and Neumann boundary conditions. A simple 
example when the term with goo in Eq. (3.49) gives contribution and the formula used in the 
papers cited above is not applicable, corresponds to an infinite plane boundary moving with 
uniform proper acceleration in the Minkowski spacetime. Assuming that the boundary is 
situated in the right Rindler wedge, in the accelerated frame it is convenient to introduce Rindler 
coordinates (r, ^,x) related to the Minkowski ones {t,x^,:x.) by transformations 



t = ^ sinh T, 



cosh r. 



(3.50) 
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and X = {x^, ... ,x^) denotes the set of coordinates parallel to the boundary. In these coordi- 
nates the metric is static, admitting the Killing vector field d/dr and has the form 



9ik = diag(^^,-l,...,-l). 



(3.51) 



The boundary dMg is defined by = a. For the boundary normal and nonzero components of 
the extrinsic curvature tensor one has 



n 



(3.52) 



Now for the components of the tensor t!^ in expression (3.32) for the surface energy-momentum 
tensor we find 



C 2 

a 



6^6l + {AC-l){l + msa/OS': 



i,k = 0,2... ,D, 



(3.53) 



and tI = 0. As we see from here, now for non-minimal coupling the surface energy density 
does not vanish for Neumann boundary condition (realized for = —(/a). An application of 
formula (3.53) will be given in Sec. 5. 



4 Integral conservation laws 

Integral conservation laws are associated with the symmetries of the background manifold. These 
symmetries are described by the Killing vector rji satisfying the equation 

ViT/fc + VkVi = 0. (4.1) 

By virtue of the symmetry of the tensor T^^°^\ from equation (3.41) one can obtain 

l=d, (v¥r(™')^Sfc) = 0. (4.2) 



Integrating this equation over the spacetime region M and using the Stoke's theorem one obtains 
the relation 



where 



JdM^ 

p(-0 = ^d^^^^,^,T(-^)«'= (4.4) 

is the volume part of the energy-momentum, and is the future pointing normal vector to 
the spatial hypersurface S. In Eq. (4.4), 7 is the determinant for the metric 7jfc induced on 
the surfaces E by the spacetime metric: 7^^ = gi\^ — UiUk- As we see from (4.3), the volume 
energy-momentum, in general, is not conserved separately. 

Projecting the Killing vector onto the boundary we find the corresponding vector for the 
boundary geometry. We will denote this vector by the index b: 

ri? = h'im- (4.5) 
The equation satisfied by this vector field is obtained as the projection of Eq. (4.1): 

+ Dkvf^ = 0. (4.6) 
In order to construct global quantities let us write equation (3.42) in the form 

DkT''' = W^mT^^"'^'^^^ . (4.7) 
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Multiplying this equation by ry|^^ , using the symmetry properties of the surface energy-momentum 
tensor and Eq. (4.6), and noting that from the orthogonality relation riir*'^ = one has 
^W^ife ^ riiT^'^, we find 

dk {VWlmr'') = V\h\hU?niT(^°'^>'^ . (4.8) 

We integrate this relation over the hypersurface dMg. After using the Stoke's theorem for the 
term on the left and noting that hl.rjf''' = rj^^ , one receives 

- = I d^x^\nd^^T^^°'^'^ . (4.9) 
JdMs 

Here we have introduced the surface energy-momentum 

p(-rf) = / d^-lx^%^,r^^ (4.10) 



9S 

where Ui is the forward pointing unit normal to the boundary in the hypersurface dMg and 
has been introduced in Sec. 3, a is the determinant of the metric (2.8) induced on 9E. 
Now summing the relations (4.3) and (4.9), we arrive at the formula 

Pe, - Psi = / d^x^n,nitnS/r(-')^*^ . (4.11) 

JdMs 



Ps = p(7°^)+p(!."'-f). (4.12) 



Here Ps is the total energy-momentum on a hypersurface E being the sum of volume and surface 
parts: 

D d(vo1) , L 

As it follows from (4.11), for any Killing vector tangent to the hypersurface dMg, n^rji = 0, the 
total energy-momentum Pj] is a conserved quantity. The subintegrand on the right of formula 
(4.11) can be interpreted in terms of the normal force acting on the boundary. This force is 
determined by the volume part of the energy-momentum tensor. The total energy-momentum 
can be also obtained from the total energy-momentum tensor (3.31) by integrating over the 
hypersurface S: 

Ps = ^ d''x^\uirikr^ . (4.13) 

To see this, first note that using the relations (see, for instance, [8]) ^/\^ = iV-^/py[, x/pyf = 
N\\/\a\ with the lapse function N and taking infinitesimaly close Si and S2, from Eq. (3.33) 
we obtain 

[ d^xy/\^\6{x;dMs)= [ d^-^x^/\^\X. (4.14) 

JT, Jd-E 

Prom the orthogonality relation njr*'^ = and the second formula in (2.1) we have 

Xuy'' = UiT^^ . (4.15) 

Hence, by taking into account (3.32) and (4.15), the integral with the surface energy-momentum 
tensor gives 

f d'^x^MuirikT^"'"'^'^'^ = I d'^-^x^\\uirik'r^^ = P^-^'^\ (4.16) 

JT, JdT 

with P^^^^^ defined as in Eq. (4.10). 

To see how the quantities P^™''* and -Pg^^^'* are related to the corresponding proper densities 
(for a similar discussion in general relativity see [6]), let us define the proper energy volume den- 
sity £^^°^\ proper momentum volume density and spatial stress p(voi)«'= on the hypersurface 
S as 
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The following relation can be easily checked: 

and, hence, expression (4.4) for the volume energy-momentum can be written in terms of proper 
quantities as 

P^°'^ = f d^x^r?ik(£(^°i)u^+p(^°')'=) . (4.19) 

Similarly, we can define the proper energy surface density e^**"''^), proper momentum surface 
density p(surf)i^ g^^^j spatial stress p(su''f)jfe ^j^g surface dH: 

with the relation 

^ ^(surf) ~i ^ p(surf)i _ (4_2l) 

The latter allows us to write expression (4.10) for the surface energy-momentum in the form 

P^^'''^=f d^-ixv1^7?fcf£(-rf)«'=+p(-'^^)*^) . (4.22) 
Jdi: ^ ' 

If a Killing vector field r^j = 77^^^^ exists which is timelike and has unit length, r}^^)^j{p = 1, 
we can take a foliation E with = r/^^^'. In this case 7y^°)p{voi)A; _ ^-^^ from (4.19) one finds 



^(vol) 



jj''x^\e^''^'\ (4.23) 



which is the total proper energy on the hypersurface S. Similarly, noting that Ukcj'i = and, 
hence, UkP^^^^^^^ = 0, from (4.22) we receive 



4-^) = / d^-^xv^^= , (4.24) 
Jdi: vl — 



where we have introduced the proper radial velocity v by using the relation ^/l + rf = 1/ VT^-^ 
(see Sec. 2). For a static boundary dT, one has n^r]^'^^ = and 7] = v = and, as it follows from 
(4.11), the total energy does not depend on the hypersurface S and is a conserved quantity. 
Note that, in general, this is not the case for volume and surface parts separately. 



5 Vacuum expectation values for the energy-momentum tensor 
on manifolds with boundaries and the vacuum energy 

In this section we will consider a quantum scalar field i^ix) on a static manifold M with boundary 
dMg. In a static coordinate frame, go/3 = and the components goo, g/sp, P, p = 1, . . . , D, do not 
depend on the time coordinate and n'' = on the boundary dMs- The quantization procedure 
more or less parallels that in Minkowski spacctimc. Let {(/?q(,t), 1^* (,t)} be a complete set of 
positive and negative frequency solutions to the field equation, obeying the boundary condition 
(3.39) on the boundary dMg. Here a corresponds to a set of quantum numbers, specifying the 
eigenfunctions. These eigenfunctions are orthonormalized in accordance with [43] 

{Va,(Pa') = -i J^d'^X^/\J\u'' [(pl{x)dk<fa' (x) - (x)9fe<y9* (x)] = Saa' , (5-1) 
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where the integration goes over a spatial hypersurface S of the manifold M with the future 
pointing normal u'^. It can be seen that scalar product (5.1) does not depend on the choice of 
the hypersurface E. Indeed, for the corresponding difference for two hypersurfaces Si and S2, 
from the field equation, by using the Stoke's theorem, one finds 

iiPa,^a'hi - (<yfa,</^a')s2 = d^Xy/\h\n^ [ip*a{x)dkiPa' (x) - (Pa' {x)dk<f*aix)] ■ (5.2) 

JdMs 

Making use boundary condition (3.39) for the functions ipa {x), we conclude that the integral on 
the right vanishes. Due to the staticity of the problem the time dependence of the eigenfunctions 
can be taken in the form e"^"*, with lo^ being the corresponding eigenfrequencies. By expanding 
the field operator over the eigenfunctions ifa{x), using the standard commutation rules and 
the definition of the vacuum state, for the vacuum expectation values of the energy-momentum 
tensor one obtains 

(0|r^,(a;)|0) = ^ T^,{^a{x),^*a{x)}, (5.3) 

a 

where |0) is the amplitude for the corresponding vacuum state, and the bilinear form Tn^{ip,ij;} 
on the right is determined by the classical energy-momentum tensor. The mode-sum on the 
right of formula (5.3) is divergent. Various regularization procedures can be used to make it 
finite (cutoff function, zeta function regularization, dimensional regularization, point-splitting 
technique). In the formulae below, it is implicitly assumed that some regularization procedure 
is used and the corresponding relations take place for regularized quantities. 

To find the total vacuum energy on M, we need the vacuum expectation values of the energy 
density (0|Tq |0). For the points x ^ dMg the latter can be evaluated by the formula 

(0|ro(^°')>) = J2[^l + {C- mrng'^^ViVm + Cr^oV/j - CRoo] (5.4) 

a 

In the static case we have 

Too = It^^'V.^oo, i?oo = --Lv^V^V^^, P,p = 1,...,D, (5.5) 
^ V^oo 



where V/3 is the covariant derivative operator associated with the Z)-dimensional spatial metric 

7^p = —gpp- For the nonzero components of the timelike Killing vector r/^°^* and the normal 
vector one has jy^*^)'^ = = I/^^qo- Substituting these into formula (4.4) and by taking 

into account relations (5.5), for the total volume energy one finds 

= E {t + is d''-^x^\nf [(C - 1/4)5^3 - C(5/3ln V5^)] <^a(x)^;(x)| , (5.6) 

where t)S is the boundary of S. As we see, the total volume energy, in general, does not 
coincide with the vacuum energy evaluated as the sum of zero-point energies for each elementary 
oscillator. The total surface energy can be found integrating the corresponding energy density 
from Eq. (3.49): 

= -E/ d''-^x^\nP[{C,-l/^)dp-adp^^Vm)]Mx)v*a{x)- (5.7) 
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For the total vacuum energy one obtains 



E = + = 1 J] Wa. (5.8) 

a 

As we have seen, the surface energy contribution is crucial to obtain the equality between the 
vacuum energy evaluated as the sum of zero-point energies of each elementary oscillator, and 
energy obtained by the integration of the corresponding energy density. 

As an example let us consider vacuum expectation value of the surface energy-momentum 
tensor for a massless scalar field induced by an infinite plane boundary moving with uniform 
proper acceleration through the Fulling- Rindler vacuum. The vacuum expectation value 
of the volume part of the energy-momentum tensor for this problem is investigated in Refs. 
[44, 45] (see Ref. [46] for the case of two parallel plates geometry). The corresponding total 
Casimir energies for a single and two plates cases are recently evaluated in Ref. [47] by using 
the zeta function technique. In the consideration below we will assume that the plate is situated 
in the right Rindler wedge and will consider the region on the right from the plate (RR region 
in terminology of Ref. [45]) corresponding to ^ > a in the Rindler coordinates (see Sec. 3). The 
boundary condition (3.39) now takes the form 

d \ 

A + a—\ip = Q, A = 2{C + msa), ^ = a. (5.9) 

For the region ^ > a a complete set of solutions to the field equation that are of positive frequency 
with respect to d/dr and bounded as ^ — >^ oo is 

ifUx) = CKUKy''^-"^\ a = {lo, k), (5.10) 

where = |k| and Kii^{x) is the MacDonald function with the imaginary order. From boundary 
condition (5.9) we find that the possible values for uj have to be roots to the equation 

AKi^{ka) + kaK[^{ka) = 0, (5.11) 

where the prime denotes the differentiation with respect to the argument. This equation has 
infinite number of real zeros. We will denote them by a; = = ujn{k), n = 1,2,.... The 
coefficient C in Eq. (5.10) is determined by the normalization condition and is equal to 

where Iioj{x) is the Bessel modified function and for a given function f{z) we use the notation 

f{z)^Af{z) + zf'{z) (5.13) 

(not to be confused with the bared notations in the conformally transformed frame in Sec. 3). 

Substituting the eigenfunctions (5.10) into the mode-sum formula (5.3) with the surface 
energy-momentum tensor and integrating over the angular part of k, for the corresponding 
vacuum expectation value one finds 



(4vr)V?(^) 

r dxx^-'f^Wrk^^ Z, A; = 0,2,..., A (5.14) 
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with r(a;) being the gamma function and (0|Tj^|0) = 0. Expression (5.14) is divergent and needs 
some regularization with the subsequent renormaUzation. One can introduce a cutoff function 
and then can apply to the sum over n the summation formula derived in Ref. [45] from the 
generalized Abel-Plana formula. An alternative approach is to insert in Eq. (5.14) a factor <^^*, 
to evaluate the resulting expression in the domain of the complex s plane where it is convergent, 
and then the result is analytically continued to the value s = 0. This approach corresponds 
to the zeta function regularization procedure (see, for instance, [48, 49] and references therein). 
The analytic continuation over s can be done by the method used in Ref. [47] for the evaluation 
of total Casimir energy. The corresponding results will be reported elsewhere. Note that on the 
plate ^ = a, energy-momentum tensor (5.14) is of perfect fluid type with the equation of state 



.(surf) 



1 + 



c 



(4C-l)(C + m.a) 



-p(surf) ^ 



(5.15) 



and effective pressure defined by (0|rf |0) = - jfpC^'irf), l,k = 2,...,D. For a minimally coupled 
scalar field this corresponds to a cosmological constant induced on the plate. 



6 Surface energy on a spherical shell 

As an application of the general prescription described in previous section, let us consider a 
spherical surface with radius a as the boundary dT, and the interior region as the manifold M in 
the case D = 3. As the vacuum expectation values for the volume part of the energy-momentum 
tensor and the total Casimir energy in this problem are well-investigated in literature (see, for 
instance, Refs. [28, 50, 51, 35]), here we will concentarte on the surface energy. Now in boundary 
condition (3.39) one has K = — 2/a, = —6[ and, hence, 

d \ 

As + a— ] (f = 0, As = AC- 2msa, r = a. (6.1) 
or J 

Inside the sphere the complete set of solutions to the field equation, regular at the origin, has 
the form 

^^{x) = ^J,{Xr)YUe,^)e-''^\ -l<n<l, / = 0,1,2,..., ^^ = 1 + ^, (6.2) 

where A = \/ uP' — m^, Jyi^z) is the Bessel function, and Yij{9,(f)) are the spherical harmonics. 
The coefficient Pa can be found from the normalization condition and is equal to 

Pa ~ ' '^'^i^) ~ T~o 9\ t9/ \ I 9 T'9/ \ ' (6.3) 

Loa [z^ — v^jJ^^z) + z^Jjfyz) 

Prom boundary condition (6.1) for eigenfunctions (6.2) one can see that the possible values for 
A have to be solutions to the following equation 

AJ^{z) + zJl{z) = Q, z = Xa, A = As-^. (6.4) 

For real v > —1 all roots of this equation are simple and real, except the case A < —p when there 
are two purely imaginary zeros. In the following consideration we will assume that A > —1/2 
with all zeros being real. Let us denote by Aj^^fe, = 1, 2, . . ., the zeros of the function AJy{z) + 
zJl{z) in the right half plane, arranged in ascending order. The corresponding eigenfrequencies 

are defined by the relation cji^^k = Note that for the massless case the existence 
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of purely imaginary zeros for Eq. (6.4) would mean that the corresponding quantum state is 
unstable. 

Substituting Eq. (6.2) into Eq. (5.7) and using the addition formula for the spherical 
harmonics, for the surface energy one obtains 

^ oo oo \2 

As it stands, the right hand side of this equation clearly diverges and needs some regularization. 
We regularize it by defining the function 

oo 

Zinis) = J^C{m)u{s), (6.6) 

where the related partial zeta function is introduced: 

k=l 

(Here and below wc use the subscripts in and cx to distinguish the quantities in the interior and 
exterior regions.) In accordance with the zeta function regularization procedure, the surface 
energy is obtained by the analytic continuation of Z\n{s) to the value s = 1/2: 

Et'^=^-^A,ZUs)l^y,. (6.8) 

The starting point of our consideration is the representation of the partial zeta function (6.7) in 
an integral form. This is done in Appendix A by making use the generalized Abel-Plana formula 
(for applications of the generalized Abel-Plana formula see also Refs. [34, 35, 36, 45, 46]). The 
corresponding integral representation for a general case of a massive scalar field is given by 
formula (A. 6). Below we will concentrate on the massless case for which from (A. 6) one has 

This integral representation is valid for 1/2 < Re s < 1. For the analytic continuation of the 
integral on the right to s = 1/2 we employ the uniform asymptotic expansions of the modified 
Bessel functions for large values of the order (see, for instance, Ref. [53]). Prom these expansions 
one has 

t ^ Wk{t) I 

where the coefficients Wj^ (t) are combinations of the corresponding coefficients in the expansions 
for I^{ux) and Il{ux). In our consideration the first three coefficients will be sufficient: 

WQ(t) = l, wi{t) = (]--A\t-]-t^, 



W2{t) =i^j-A + A'jt'+^A-:^jt^ + -^^ (6.11) 
Note that the functions Wk{t) have the structure 

k 
j=0 
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where the coefficients w^j depend only on A. Following the usual procedure applied in the 
analogous calculations for the total Casimir energy (see, for instance, Refs. [50, 52, 28]), we 
subtract and add to the integrand in Eq. (6.9) N leading terms of the corresponding asymptotic 
expansion and exactly integrate the asymptotic part. By this way Eq. (6.9) may be split into 
the following pieces 

(6.13) 



where 



C(in).(^) = drl(^) + Cv(^)' 



(1) 



■'{m)u 



''(in)i^V*; 



- .2(1-.) £^ 



ivrs 

TT Jo 



dxx 



l-2s 



N 

E 

fc=0 



tWk{t) 



2(l-s) 



sinvrs 



TT 



dxx 



l-2s 



luii^x) 



N 



-E 



,,fc+i 



fe=0 



(6.14) 



(6.15) 



The number of terms N in these formulae is determined by the condition of convergence for the 

(s) at s = 1/2. It can be seen that it is sufficient to take N > 2. For larger 



expression 2^;^q vQ 



N the convergence is stronger. Now by taking into account Eq. (6.12) and the expression for t, 
the integral in Eq. (6.14) is easily evaluated in terms of gamma function. For the corresponding 
sum over I this gives: 



1 



1=0 



N 

2f(i) 



Ch 2s + A; - 2, 



1 



V- r(j + 

j=0 



ri + .) 



J + 



Mi) 



(6.16) 



where we have introduced Hurwitz zcta function Ch(S)O)- In the following calculations we will 
take the minimal value N = 2. The expression on the right of Eq. (6.16) has a simple pole at 
s = 1/2 coming from the poles of Hurwitz and gamma functions. Laurent-expanding near this 
pole one receives 



where 



,(as) 



is) 



Z 



(as) 



s - 1/2 



4^'^ + 0{s - 1/2), 



(6.17) 



7(as) 



^(as) 



1 

2^ 
1 

TT 



A' 



'-A + 1- 
3 20 



A 



A^ (7 + 31n2) + - (l-7-31n2) 



2^^^ ' 120 



- (0.0382 - 0.552A + 2.657^^) , 

TT ^ ^ 



(6.18) 



(6.19) 



with the Euler's constant gamma and Riemann zeta function Cr('S)- 

As regards to Eq. (6.15), for A/" > 2 it is finite at s = 1/2 (the series over I converges as 
1//^) and we can directly put this value and evaluate the integral and sum numerically. Now 
the expression for the surface energy contains pole and finite parts: 



p(surf) _ 77i(surf) p(surf) 
^in - ^{m)p -^(in)f ' 



where for the separate contributions one has 



^(surf) 
(in)p 



1-4C ^^if 
a s-1/2' 



(m)f 



1-4C 



A, 



1=0 



(6.20) 



(6.21) 
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It is well-known that the total vacuum energy inside a spherical shell also contains pole and finite 
contributions [28, 50, 51]. In particular, the asymptotic part of the total Casimir energy inside 
a spherical shell has a structure similar to Eq. (6.16) and the corresponding pole part comes 
again from the poles of Hurwitz and gamma functions. Note that taking interior and exterior 
space together, in odd spatial dimensions the pole parts cancel and a finite result emerges. Now 
let us see that the similar situation takes place for the surface energy. 

Consider the scalar vacuum outside a spherical shell on which the field satisfies boundary 
condition (6.1). To deal with discrete spectrum, we can introduce a larger sphere with radius 
R, enclosing one of radius a, on whose surface we impose boundary conditions as well. After 
the construction of the corresponding partial zeta function we take the limit i? — ^ oo (for details 
of this procedure see Ref. [50] in the case of the exterior Casimir energy and Ref. [35] for the 
corresponding vacuum densities). As a result for the surface energy in the exterior region one 
obtains 



Et''^=^-^AsZUs)\s=V2. (6.22) 



with 



ZUS) = E^C(ex).(^), C(ex).(^) = ^ rdxx'-^^^ (6.23) 

In deriving Eq. (6.23) we have assumed the condition ^ < 1/2 under which the subintegrand 

has no poles on the positive real axis. The corresponding expression differs from the interior 
one by the replacement ly Ky. This relation between the interior and exterior quantities is 
well-known for the total Casimir energy and vacuum densities. Now the analytical continuation 
needed can be done by the way described above for the interior case using the uniform asymptotic 
expansions for the function Ky{x). The corresponding formulae will differ from Eqs. (6.14)— 
(6.16) by replacements —>■ and Wk{t) —>■ {—l)''~^^Wk{t). Now we see that in the sum 

z(^){s) = Z^^\s) + zi^Hs), (6.24) 

the summands with even k and, hence, the poles terms cancel (recall that the pole terms in 
the separate interior and exterior surface energies are contained in k = and k = 2 terms). 
Moreover, for s = 1/2 the term with A; = 1 in the expression for Z^^\s) vanishes as well due 
to Ch(0, 1/2) = 0. In particular, the asymptotic part of the total surface energy vanishes for 
N = 2. Taking this value and using the expression of the function wi{t), for the total surface 
energy one obtains 



p;(surf) ^ ^(surf) E;(surf) 
-^in ~ ex 

Jo [lyix) Ky{x) x^ + v^\x^ + v^ 



na 

1=0 



, (6.25) 



where we have changed the integration variable ux ^ x. In Fig. 2 we have plotted the quantity 
as a function on the Robin coefficient As- This function vanishes at As ~ 0.37, 
has a maximum at Ag ~ 0.07 and a local minimum at As ~ 0.54. 



7 Conclusion 

On background of manifolds with boundaries the physical quantities, in general, will receive 
both volume and surface contributions and the surface terms play an important role in various 
branches of physics. In particular, surface counterterms introduced to renormalize the diver- 
gencies in the quasilocal definitions of the energy for the gravitational field and in quantum 
field theory with boundaries are of particular interest. In the present paper, continuing the 
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Figure 2: The total surface energy for a spherical shell, aE(^"''f)/(l - 4C), as a function on the 
coefficient As in boundary condition (6.1). 

investigations started in Refs. [34, 35, 36], we argued that already at classical level there is a 
contribution to the energy-momentum tensor of a scalar field, located on the boundary. We 
have considered general bulk and boundary geometries and an arbitrary value of the curvature 
coupling parameter for the scalar field. Our starting point is the action (2.2) with bulk and 
surface contributions (2.3) and (2.4) (see also (2.7) for the corner term). The surface action con- 
tains the generalization of the standard Gibbons-Hawking surface term for the problem under 
consideration, and an additional quadratic term. We evaluate the energy-momentum tensor in 
accordance with the standard procedure, as the functional derivative of the action with respect 
to the metric. In the variational procedure we held fixed the metric on initial and final spacelike 
hypersurfaces Si and 112 (see Fig. 1) and allowed variations on the timelike boundary dMg. 
The surface energy-momentum tensor is obtained as the functional derivative with respect to 
the variations of the metric tensor on SM^. As the various surface terms arising in the variation 
procedure are essential for our consideration, the corresponding calculations are presented in 
detail for the bulk and boundary parts of the action in Sec. 3. The bulk part of the energy- 
momentum tensor has the well-known form (3.11) and the surface energy-momentum tensor is 
given by expression (3.32) with Tik from Eq. (3.30), and the delta function locates this part 
on the boundary dMg. The surface part of the energy- momentum tensor satisfies orthogonality 
relation (3.34) and does not depend on the foliation of spacetime into spacelike hypersurfaces. 

Similar to the metric variation, taking the variation of the action induced by an infinitesimal 
variation in the field, we have fixed the field on the initial and final spacelike hypersurfaces Ei 
and T,2 and allowed variations on the hypersurface dMs- From the extremality condition for the 
action under the variations of the field on dMg we have received the boundary condition on dMg, 
given by Eq. (3.39). On the class of solutions to the field equation, the volume part of the energy- 
momentum tensor satisfies the standard covariant conservation equation (3.41), and for the 
surface part we have obtained equation (3.42) which describes the covariant conservation of the 
surface energy-momentum tensor up to the flow of the energy-momentum across the boundary 
into the bulk. For a conformally coupled massless scalar field both volume and surface energy- 
momentum tensors are traceless and the boundary condition for the field on dMg coincides 
with the conformally invariant Hawking boundary condition. Further we have considered the 
surface energy density for a static manifold with static boundary and have shown that the 
expression for this quantity, previously used in Refs. [40, 34, 35, 36], is obtained as a special 
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case with ^fQO = 1- As an example of geometry when the formula used in the papers cited 
above needs a generalization, we considered a plate moving with uniform proper acceleration 
in the Minkowski bulk. In Sec. 4 the integral conservation laws are investigated which follow 
from the covariant differential conservation equations (3.41) and (3.42). Under the presence 
of Killing vectors the corresponding volume and surface conserved charges are constructed by 
formulae (4.4) and (4.10). The total energy momentum is obtained as a sum of these parts. The 
difference of the total energy-momentum for two spacelike hypersurfaces Si and S2 is related 
to the projected onto normal part of the volume energy-momentum tensor by Eq. (4.11). 
The subintegrand on the right of this formula can be interpreted in terms of the normal force 
acting on the constraining boundary. This force is determined by the volume part of the energy- 
momentum tensor only. Further, we have shown that the total energy-momentum can be also 
obtained by integrating the total energy- momentum tensor over a spacelike hypersurface S. We 
have expressed the volume and surface energy-momenta in terms of the proper energy volume 
and surface densities and proper momentum volume and surface densities, formulae (4.19) and 
(4.22). A special case is considered when a timelike Killing vector exists for the bulk metric. 
By choosing an appropriate foliation we have seen that in this case the corresponding volume 
energy-momentum coincides with the integrated volume proper energy. For the general case of 
the non-orthogonal foliation the integrated surface energy is given by formula (4.24). In the 
case of an orthogonal foliation, which corresponds to a static boundary, the total energy is a 
conserved quantity. We have emphasized that, in general, this is not the case for separate volume 
and surface parts. 

In Sec. 5 we have considered a massive quantum scalar field on a static manifold with static 
boundary. By integrating the vacuum expectation values of the energy densities we have found 
the volume and surface parts of the vacuum energy. We have shown that (i) the integrated 
volume energy, in general, differs from the vacuum energy evaluated as the sum of zero-point 
energies for each normal mode and (ii) this difference is due to the additional contribution coming 
from the surface part of the vacuum energy. These results for the special cases of flat, spherical 
and cylindrical Robin boundaries on the Minkowski bulk have been discussed previously in 
Refs. [34, 35, 36]. As an application of general formulae we have derived a formal expression 
for the surface energy-momentum tensor induced on an infinite plane boundary moving with 
uniform proper acceleration through the Fulling-Rindler vacuum. This tensor is of barotropic 
perfect fluid type with the equation of state (5.15) and in the case of a minimally coupled scalar 
corresponds to a cosmological constant on the plate. 

As an application of a general formula for the surface energy-momentum tensor, in Sec. 6 
we have evaluated the vacuum expectation value of the surface energy for a quantum scalar 
field for interior and exterior regions of a spherical shell. As a regularization procedure we have 
used the zeta function technique. An integral representation for the related partial zeta function 
well suited for the analytic continuation is constructed in Appendix A, by using the generalized 
Abel-Plana summation formula over zeros of a combination of the Bessel function with deriva- 
tives. Subtracting and adding to the integrands the leading terms of the corresponding uniform 
asymptotic expansions, we present the interior and exterior surface energies as sums of two 
parts. The first ones are finite at s = 1/2 and can be easily evaluated numerically. In the sec- 
ond, asymptotic, parts the pole contributions are given explicitcly in terms Hurwitz and gamma 
functions. The remained pole terms are characteristic feature of the zeta function regularization 
method and have been found for many cases of boundary geometries. For an infinitely thin shell 
taking interior and exterior regions together, the pole parts cancel and the total surface energy 
is finite and can be directly evaluated by making use formula (6.25). The cancellation of the pole 
terms coming from oppositely oriented faces of infinitely thin smooth boundaries takes place in 
very many situations encountered in the literature. It is a simple consequence of the fact that 
the second fundamental forms are equal and opposite on the two faces of each boundary and. 
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consequently, the value of the corresponding coefficient in the heat kernel expansion summed 
over the two faces of each boundary vanishes [48] . In even spatial dimensions there is no such a 
cancellation. 

In the consideration above we have assumed a smooth timelike boundary dMg- For non- 
smooth boundaries the corner terms in the action will be important. In this case, in the variation 
of the action with respect to the metric, terms will appear with the integration over the corners. 
This terms will lead to the additional part in the energy-momentum tensor located on the 
corners. 
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A Integral representation for the partial zeta function 



In this appendix we will derive an integral representation of partial zeta function (6.7). For this 
we use the summation formula over zeros Xy^k derived from the generalized Abel-Plana formula 
[54, 55] (see Eq. (3.26) in Ref. [55]): 



k=l 





1 



dxf(x)Jl{x) 
KAx) 



^ I dx--^ll{x) /(a;e-/2) + /(a;e— /2) , (A.l) 



where the function Ti,{z) is defined by Eq. (6.3) and we use the notation indroduced in Eq. 
(5.13). This formula is valid for functions f{z) analytic in the right half plane and satisfying 
the condition 

!/(z)|<^, a>0, l^l^oo, (A.2) 



with a constant AIq. To transform the first integral on the right of formula (A.l), let us apply 
this formula to the case ^ = (the sum is over zeros of the Bessel function Jy{z)). Now the left 
hand side is zero and from (A.l) one receives 



1 /"°° r 
dxf{x)Jl{x) = - / dxU{x)Ky{x) f{xe^''^) + f{xe-^''^) 

7o L 



The resubstitution into Eq. (A.l) yields 



K,kf{K,k) _ 1 r, ijAx)_ 



f{xe^^l^) + f{xe 



-7r^/2^ 



(A.3) 



(A.4) 



where we have used the relation Ty{\v,k)Jv{^i',k) = ^u,k/{^u fe + ~ Wronskian for 

the Bessel modified functions. As a function f{z) in this formula let us choose 



which satisfies the condition (A.2) for Re s > 1/2. We obtain 



^l,ki'^l,k +/^^) 



smvrs 



vr 



dx 



X 



h{x) {x^-fi^y 



(A.5) 



(A.6) 



This formula gives an integral representation of the partial zeta function introduced in section 6. 
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